Spin Gap in the Zigzag Spin-1/2 Chain Cuprate Sro.gCao.iCuC^ 
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We report a comparative study of 63 Cu Nuclear Magnetic Resonance spin lattice relaxation rates, 
Tj - , on undoped SrCu02 and Ca-doped Sro.gCao.iCuC^ spin chain compounds. A temperature 
independent T-[~ is observed for SrCu02 as expected for an S = 1/2 Heisenberg chain. Surprisingly, 
we observe an exponential decrease of T{~ for T < 90 K in the Ca-doped sample evidencing the 
opening of a spin gap. The data analysis within the J1-J2 Heisenberg model employing density- 
matrix renormalization group calculations suggests an impurity driven small alternation of the 
J2-exchange coupling as a possible cause of the spin gap. 

PACS numbers: 76.60.Es, 75.40.Gb, 75.10.Pq, 75.40.Mg 



One-dimensional (ID) quantum magnets are intrigu- 
ing model systems for studying correlated many-body 
quantum physics. They are theoretically well treat- 
able which often yields clear-cut predictions about their 
ground states and excitation spectra (see, e.g., [l|). A 
particularly important case is the exactly solvable uni- 
form antiferromagnetic (AFM) 5 = 1/2 Heisenberg chain 
model (HM) which exhibits a gapless spectrum of elemen- 
tary spinon excitations and a ground state lacking long- 
range magnetic order. Spin gaps in the excitation spectra 
of closely related quantum systems arise only under cer- 
tain circumstances. For instance they occur naturally in 
systems with a pronounced chemically caused dimerized 
structure or in spin-Peierls compounds with a sponta- 
neously dimerized ground state @, Q . The latter become 
highly nontrivial in the presence of frustrated exchange 
interactions, being therefore a hot topic in the field of 
quantum magnetism. The parent compound SrCu02 
considered here, where Cu 2+ (3c? 9 , S = 1/2) ions bridged 
via oxygen ligands form well defined chain units along the 
crystallographic c-direction (see inset of Fig. Q]), however, 
does not belong to either of them. It is commonly con- 
sidered as one of the best realizations of the ID 5 = 1/2 
AFM HM. Due to a very small and frustrated coupling 
between different zigzag chains it shows a Neel-like or- 
dering below T N ps 2K only, despite the huge effec- 
tive AFM next nearest neighbor (NNN) exchange integral 
J 2 ~ 2000 K The nearest neighbor (NN) coupling 

J\ is much weaker, ferromagnetic (FM) and frustrated 
with an estimated ratio of \Ji/J?\ < 0.1 — 0.2 [f|. In- 
elastic neutron scattering reports the onset of a weak 
sublattice magnetization below T cl = 5K, static on a 
ns time scale, but the absence of 3D longe range order 
for T > 10~ 4 J/ks- Furthermore, it suggests a gapless 
spinon excitation spectrum for this compound down to 
the lowest measured energy transfer of 0.5 meV How- 
ever, despite inevitably present spin-phonon coupling, a 
spin-Peierls transition has not been found at low temper- 
ature, in contrast to the well-known CuGeOa [Io| . 

In this Letter we address this dichotomy by studying 
the effect of isovalent Ca substitution on the Sr sites in 



SrCu02, outside the spin chain unit. Our comparative 
Nuclear Magnetic Resonance (NMR) study investigates 
the temperature dependence of the 63 Cu nuclear spin lat- 
tice relaxation rate T^ 1 for single crystals of SrCu02 and 
Sro.gCao.iCu02. While is temperature independent 
in SrCu02 in compliance with the 5=1/2 HM and the ob- 
served gapless spinon continuum an exponential de- 
crease of Tf 1 below 90 K for Sro.9Cao.iCu02 evidences 
a significant gap in the spin excitation spectrum. Wc 
present a theoretical analysis of this striking observation 
in terms of the J1-J2 ID HM for two possible scenarios: 

(i) the AFM- AFM frustrated HM for which, based on 
qualitative arguments only, a small spin gap has been pre- 
dicted in the crossover region between a commensurate 
dimerized-like and an incommensurate spiral-like state 
[ll[; (ii) the standard spin-Peierls scenario with an al- 
ternation of the largest AFM NNN coupling J2 in the 
presence of a realistic relatively small NN coupling J\ of 
arbitrary sign. We conclude on the relevance of scenario 

(ii) for the occurrence of a spin gap in Sro.gCao.iCu02 
where subtle structural distortions due to Ca substitu- 
tion might cause a small modulation of J2 and trigger a 
drastic change of the spin excitation spectrum. 

Single crystals of SrCu02 and Sro.gCao.iCu02 were 
grown by the travelling solvent floating zone method. 
Their high quality, i.e. single crystallinity, stoichiome- 
try, the absence of phase irregularities, and second phase 
inclusions, was verified by x-ray diffraction, polarized 
optical microscopy, energy-dispersive x-ray spectroscopy 
and inductively coupled plasma optical emission spec- 
trometry [T^]. The Ca-content of Sr .9Ca .iCuO2 was 
found to be (10.0 ± 0.2) %. Structural refinement on 
Sro.gCao.iCu02 showed no measurable Ca occupancy of 
the Cu sites. The averaged change of the Cu-O-Cu bond 
angle and bond distance along the direction of J2 is only 
0.35 %. The 63 Cu NMR spin lattice relaxation rate, Tf 1 , 
was measured by inversion recovery in a magnetic field of 
H = 7.0494 T for H \\ a, b and c, respectively The values 
of Tf 1 were obtained by fitting the recovery curves to the 
standard expression for magnetic relaxation of a nuclear 
spin 7 = 3/2 13|, H| yielding one single T\ component 



2 



2 3 



v 2 
I- 

1 


7 

6 

5 

v w 4 

"E 3 

V 2 
1 



■u. ■ * * i 



(a) 
SrCuO, 



- H || c 

• H || b 

+ H || a (x6 



■ i i 1 ! 



»* 4 



(b) 



Sr Ca CuO 

0.9 0.1 2 



50 100 150 200 250 300 
T(K) 



10000 



1000 



100 r 




0.14 



FIG. 2: (Color online) Spin lattice relaxation rate Tj - of 
Sro.gCao.iCu02 vs. inverse temperature. Lines are fits of 
Eq. © with a fixed spin gap of A = 50 K at low temperature 
(10 - 40 K). For clearness T^~* !f has been subtracted from the 
data and the fits. Inset: the same for H\\b in H = 2T. 



FIG. 1: (color online) Spin lattice relaxation rate T^ 1 of 
SrCu02 (a) and Sro.gCao.iCu02 (b) for different field direc- 
tions. Inset: zigzag structure and main couplings. 



at any temperature |15[. 

For SrCu0 2 , T^ 1 remains constant over the whole tem- 
perature range, and increases only slightly at low temper- 
ature [Fig. [TJa)] [l6j . Note that the error bars increase 
with decreasing temperature due to a strong broadening 
of the NMR resonance lines at low temperature 17| . The 
observed behavior strongly supports the gapless excita- 
tion spectrum in SrCu0 2 and is consistent with theoret- 
ical calculations for the 5=1/2 AFM HM predicting a 
temperature independent T^ 1 18-2(]|- Fig. [TJb) shows 
the temperature dependence of T^ 1 for Sro.gCao.iCu0 2 - 
Similar to the undoped sample, T^ 1 shows the expected 
temperature independent behavior for T > 90 K, but un- 
expectedly decreases exponentially by more than one or- 
der of magnitude for T < 90 K 

In case of magnetic relaxation the spin lattice relax- 
ation rate probes directly the imaginary part of the dy- 
namic susceptibility x" of the electronic spin system: 



T 



const 



x"(<?» 



(1) 



Here A± is the hypcrfinc coupling, q the wave vector 
and u> the NMR frequency. The observed exponential 
decrease of T^ 1 gives therefore evidence for the opening 
of a gap A in the excitation spectrum of Sro.gCao.iCu0 2 
that is manifested by an exponential reduction of x" ■ We 
estimated the magnitude of the gap by fitting our data 
to an activated temperature dependence: 

S (-A/T) 



which applies for T < A. To account for the non- 
vanishing value of T^ 1 as T — > 0, we added a constant 
offset T£" oH to the activation behavior. The results shown 
in Fig. [2] clearly demonstrate that for T < 40 K the de- 
crease of T^ 1 can be well described with an isotropic 
spin gap of A = 50 K. A possible field dependence could 
be excluded by supplementary measurements of T^ 1 in 
H = 2T, also yielding A = 50 K (see inset of Fig. 

The opening of a spin gap suggests that the magnetic 
properties of the Cu spin chains change drastically by 
substituting Sr with isovalent Ca. Note that neither in 
susceptibility nor in other thermodynamic properties any 
indication for a spin gap was found. This is not surprising 
since owing to the very large value of J2 the intrinsic 
static magnetic susceptibility at T -C J2 is very small. 
Therefore the thermodynamics might just overlook the 
gap since the total static magnetic response at low T is 
dominated by other contributions, such as the Van-Vleck 
susceptibility and Curie-like impurity contributions. 

Two different approaches to explain such a doping- 
induced spin gap in Sro.gCao.iCu0 2 will be discussed 
based on the following spin Hamiltonian written in the 
adiabatic approximation: 



H=j 1 j2s i - s l+1 + j 2 5^(1 + (-lrs)^ ■ s,, 



i + 2- 



(3) 



Tf L {T) = Tf o ^ + const 



(2) 



Here Si is a spin-i operator at site i and a 
\/2sin |(2i - 1), Ji and J 2 are the NN and the NNN ex- 
change interactions and S is the modulation of J 2 - The 
spin gap is evaluated as the energy difference between 
the lowest triplet state and the singlet ground state, 

A(L) = E^L) - E (L), A = lim A(Z), (4) 
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well fitted by the expression 



FIG. 3: Spin gap as a function of J2/J1 in the AFM-AFM sec- 
tor of the J1-J2 HM. The experimental value of the spin gap in 
units of Ji is depicted by dashed lines adopting J2 = 2000 K. 



where E n (L) is the n-th eigenenergy (n = denotes the 
ground state) of the system with length L. We employ 
the density-matrix renormalization group (DMRG) tech- 
nique, a powerful numerical method for various ID quan- 
Open-end boundary conditions are ap- 
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turn systems 

plied in the chain direction and it enables us to calculate 
the energies of the ground-state and low-lying excited- 
states quite accurately for very large but always finite 
systems (up to the order of ~ 10 3 sites). 

Scenario (i) gives attention to the frustrated NN cou- 
pling J\. For FM NN interaction Jx < 0, whose absolute 
value is much smaller than the NNN interaction J 2 , as 
it is likely the case in SrCuOo, the zigzag chain is not 
expected to show a spin gap However, in the case 
of a small AFM NN coupling J\ > 0, the zigzag chain 
should exhibit a spin gap which scales exponentially with 
the ratio of the two competing couplings [H( : 



A oc exp —const 



■h 
Ji 



(5) 



According to the Goodenough-Kanamori-Anderson 
(GKA) rules (H a Cu-O-Cu bond angle of 180° yields a 
strong AFM exchange between the Cu spins (J 2 ), while 
an angle of 90° results in a weak FM coupling (Ji). 
However, as soon as the bond angle slightly deviates 
from 90° , the interaction J\ can become AFM again [24| . 
Local structural distortions in Sro.gCao.iCu02 induced 
by Ca ions could yield such a deviation, resulting in a 
local spin gap corresponding to Eq. (|5|). 

Within this scenario (6 = and J\ > 0) wc study 
chains with several lengths L = 512 to 4096, keeping m = 
1200 density- matrix eigenstates in the renormalization 
procedure. In this way, the discarded weight, Wd, is less 
than 10 -7 , while the maximum error AE in the ground- 
state and low-lying excited states energies is less than 
10~ 7 — 10 ~ 6 . The extrapolated results of the spin gap 
are shown in Fig. [3] For J 2 /Ji > 2, the spin gap can be 



A/Ji = 1.6 exp (-1.7 JaM), 



(6) 



which might be useful for the description of frustrated 
systems in general. However, for the zigzag chain under 
consideration with the experimentally obtained gap of 
A = 50 K a value of J\ ~ 1000 K would be required 
to explain our NMR data. This is unrealistically large 
for a Cu-0 bond angle only slightly below or above 90° 



24, 



25| . Hence, scenario (i) can be discarded. 
Scenario (ii) examines the possibility of an alternating 
AFM NNN exchange J 2 (l ± S), which is known to cause 
an intrinsic spin gap. Within this scenario we assume 
that the Ca substitution might induce either directly or 
via a softening of the lattice and a concomitant strongly 
increased spin-phonon interaction a modulation of the 
NNN coupling constant J 2 , whose simplest realization is 
a periodic alternation. The validity of this scenario will 
be examined for FM and AFM Ji. 

We study chains with several lengths L = 64 to 392, 
keeping m = 2000, which ensures Wd < 10~ 14 and AE < 
10~ 12 . In Fig.HJa), the spin gap is plotted as a function 
of 5 for Ji = 0. We estimate S = 0.0027 to reproduce the 
experimental value of the spin gap A = 50 K for J 2 = 
2000 K. Next, we introduce the NN exchange interaction 
J\. The spin gap is calculated as a function of J\ with 
fixed 6 = 0.0027 [Fig. gfb)]. We find that the value of A 
is hardly affected by J\. Thus, only a small alternation 
factor 5 = 0.0027 is needed to reproduce the spin gap 
found by NMR, irrespective of the strength and sign of 
the interchain coupling \J± \ <C J 2 . 

Beyond the simple alternation of J 2 , a more realistic 
situation would be given by including also the direct im- 
pact of disorder in the exchange coupling caused by the 
Ca impurities. Thus, we consider an additional term to 
Eq. (O, Hdi S = 2D YliJi £ iSi ■ Si+2, where is defined 
by a box probability distribution V(ei) = 0(1/2 — \Ei\) 
with the step function 0(x) and the disorder strength is 
controlled by D. We calculate the spin structure factor 
S(q = 7r, oj) on random sampling 30 realizations of V(ei) 
for L — 384, using dynamical DMRG [26|, and take an 
average of the results for each D. Fig. [4jc) shows the 
averaged S(q = 7r,u>) for several D values with S = 0.1 
and Ji = 0. We find that the delta peak for D = is 
broadened by the disorder and a "pseudo spin gap" -like 
behavior appears. Note that a relatively large 8 = 0.1 
is chosen to be able to clearly follow the D-dependence. 
For D > S the singlet pairs are partially collapsed, i.e., 
magnetic moments are locally induced, and the spectral 
weight at u> = becomes finite. This must be related 
to the occurrence of a finite, doping-dependent residual 
relaxation rate T-f^ [see Eq. ([2]) and Fig.QJb)]. 

In summary, measurements of the 63 Cu NMR spin lat- 
tice relaxation rate Tf 1 in a single crystal of the S — 1/2 
chain compound Sro.gCao.iCu0 2 reveal a striking expo- 
nential decay of T^ 1 below 90 K in sharp contrast to a 
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FIG. 4: The same as in Fig. [3] for an alternating AFM NNN 
exchange J2(l ± 5) at (a) zero and (b) weak arbitrary NN- 
exchange J\. (c) Spin structure factor for several disorder 
strengths D with 8 = 0.1 and Ji = 0. The inset of (a) shows 
(A/J2) 1 ' 5 for small alterations 8. 

constant T^ 1 in the parent compound SrCuC>2 [13] • The 
observed decay points to the opening of a spin gap with 
A = 50 K. Our DMRG calculations show that this spin 
gap can be well reproduced within the J\-Ji ID HM con- 
sidering an alternating AFM NNN J^il ± 6) with only a 
small alternation of about 5 = 0.0027, regardless of the 
strength and sign of the NN J\ . Such a remarkable sen- 
sitivity of the low energy excitation spectrum of the S = 
1/2 chains in SrCu02 to structural substitutional defects 
outside the zigzag chain appears to be crucial for the 
understanding of the spin dynamics in low dimensional 
cuprates. Our results call therefore for extensive experi- 
mental and theoretical work on this and related systems 
to understand the impact of additional degrees of free- 
dom on the fundamental magnetic phenomena occurring 
in quantum spin magnets on the basis of cuprates. For 
example diffraction experiments may confirm the dimer- 
ization j29| . while neutron scattering, resonant inelastic 
x-ray spectroscopy, optics and specific heat may observe 
the spin gap and/or a possible phonon softening. So far 
there seem to be no generally accepted microscopic mod- 
els for spin-Pcicrls physics in cuprates. We hope that our 
work will stimulate the development of such models. 
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